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ji^. Abstract 

. A reduction of the Dirac-Maxwell equations in the case of static cylindrical 

symmetry is performed. The behaviour of the resulting system of o.d.e.s is examined 
^ ' analytically and numerical solutions presented. There are two classes of solutions. 



The first type of solution is a Dirac field surrounding a charged "wire." The 
Dirac field is highly localised, being concentrated in cylindrical shells about 
the wire. A comparison with the usual linearized theory demonstrates that 
this localisation is entirely due to the non-linearities in the equations which 
result from the inclusion of the "self-field" . 

The second class of solutions have the electrostatic potential finite along the 
axis of symmetry but unbounded at large distances from the axis. 
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1 Introduction 



This paper is concerned with the system of partial differential equations known as the 
Dirac-Maxwell equations; the Dirac field describing spin-half matter and the Maxwell field 
mediating the electromagnetic interactions of the Dirac field. In an earlier paper, one 
of us (C.J.R.) used the two-spinor formalism to rewrite the equations in a novel way. In 
fact the electromagnetic potential can be eliminated from the equations altogether. (For 
further background into the two-spinor formulation of the Dirac equation, see 

In a reduction of the equations was performed in the static case (see section 2, 
below) and a signifigantly simplified set of (nonlinear) p.d.e.s presented. Finally, after 
imposing spherical symmetry, the equations were solved (numerically, at least). Unusual 
features to emerge from that investigation were the existence of highly compact objects 
and a magnetic monopole. In the present paper we will solve (numerically) the cylindri- 
cally symmetric, static Dirac-Maxwell equations. In this case we find no monopole-like 
solutions. Perhaps this is not unexpected, since there is slightly more regularity in a 
cylindrical system as opposed to a spherical system - ln(r) behaviour as distinct from 
1/r behaviour. However, we do find highly localised solutions - cylindrical shells about 
a central charged "wire" . When we compare this to the Dirac field solved in an external 
potential of the same type, no such localisation is apparent. 

The cylindrical reduction is also distinguished from the spherical case by the existence 
of unbounded solutions - solutions in which the Maxwell potential is unbounded at large 
distances from the symmetry axis. These solutions are also unbounded in the total charge 
(per unit length of "wire" ) . 

Very few explicit solutions to the full Dirac-Maxwell equations are known and all of 
these are numerical solutions (or partly so) - see and IQ. Work on the existence 
theory has progressed steadily, with some important recent work (see and 0) providing 
local existence for soliton like solutions. 

The solutions in and those presented here show that these equations are capable 
of representing highly localised structured, entities. Such behaviour is not even hinted at 
when one examines the usual linearised theory (see section5.2), in which one ignores the 
Dirac current as a source for the Maxwell field. 
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2 The Static Dirac- Maxwell Equations 



The Dirac-Maxwell equations in standard notation are, 

'j°'{da — ieAa)ip + 'i^i^ = 0, where a = 0,...,3 (1) 

and = -47rej/3, (2) 

where = A^^^ - ^a,/3 (3) 

with the current density given by 

J a = ^lai^- (4) 

Following [cjr] we use the 75 diagonal or van der Waerden description (2-spinor notation, 
see PI). The Dirac 4-spinor (or bispinor) is 



The Dirac equations are 



v2 

{d^^ + teA^^)vA+'-^u^ = 0. (5) 
V2 

In [cjr] it was shown that these equations can be solved for the electromagnetic potential, 
provided the 2-spinors meet a non- degeneracy condition, w^va 7^ 0. The result is, 

A^^ = ( V^d^^UB + U^d^^VB + '-^{u^u^ + v^v^)\ . (6) 

e{u^vc) V V2 / 

We are interested in the classical Dirac-Maxwell equations, so we will also require that 
A^^ is a real vector field (see [cjr]). This leads to the following first order differential 
equations for ua and va ( "reality conditions" ) : 

d^'^iuAUA) = -^(mSc-m%) 

d^'^ivAVA) = ^(nSc-u%) 

UAd^^^V^-V^d^^UA = 0. (7) 

The expression for the potential (6), the reality conditions (7) and the Maxwell equa- 
tions (2), now constitute the full set of nonlinear partial differential equations for the 
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Dirac-Maxwell system. We still have, of course, the U{1) gauge freedom. 



Aa ^ Aa + \daC,-i a direct consequence of (6). 

Following [cjr] we now impose the static condition: there exists a Lorentz frame in 
which there is no current "flow" , j"' — S^j^ . Under this condition we have 

= e*^V2a°^^M^ , (8) 

where x is an arbitrary real function and a"^^ are the three Pauli matrices with the 2x2 
identity matrix as the zeroth a matrix (van der Waerden symbols). 
The gauge can be fixed by the choice 

y} = ygtM , (9) 

where 77, X and Y are real functions. 

The static Dirac-Maxwell equations can now be written down, we follow [cjr] and write 
them in 3-vector notation. This is done by first introducing vectors V and A, 



V = (2XFcosr/, 2XYsmr], X^ - 
A^{A\A^, A"^), where = a'^^A 



BB 



Our Maxwell-Dirac equations are now given as follows. 
The electromagnetic potential is 

e ^ 2e(X2 + Y^) dt 2e{X^ + Y^) 



A 



ot 



(10) 



2e(X2 + F2) 
The reality conditions are 

^{X' + Y') = (11) 
V.V = -2m{X^ + Y^)smx (12) 
^ + (Vx)xV = 0. (13) 

Together with the Maxwell equations for A with current vector j, as above. 
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3 Cylindrical Symmetry 



The Dirac field given by the 2-spinors ua and v"^ determines two gauge invariant null 
vectors P = a'^^'^UAUj^ and = a'^^'^VAVA- f^ct, 

{n = ^{x' + Y',v) and K) = i= (x^ + ^ -V ) , 

assuming that the field is static. Note that the current vector = l'^ + n°'. The first 
reality condition tells us that the zeroth component of and are independent of time. 
We also assume that and n" are independent of z in keeping with the assumption of 
cylindrical symmetry. Here (p, 0, z) are the cylindrical polar coordinates corresponding 
to our original Cartesian coordinates. 

Applying these conditions one finds that rj = (p, X = X{p) and Y = Y{p). The second 
of the reality conditions then implies that x = x{p)- The third reality condition then 
gives 

dp 

If we take x = constant then we find = constant, the Maxwell equation then implies 

that X = Y = 0. However, the non- degeneracy condition, u^va 0, excludes this 

possibility. So we assume x'(p) from now on and consequently, X = Y. Our equation 

for the electromagnetic potential ([T0|) , shows that the vector potential A vanishes. There 

is no magnetic monopole, as distinct from the spherical case. 

The equations now simplify to the following ordinary differential equations: 

m 1 
Ao = —cosx+^-r- 
e 2e dp 

A (px"^) = -2mpX^ sin x 
dp ^ ' 

± (.f) - S.V2eX', (14) 
Introducing dimensionless variables 



A = 


—Ao 




m 


P = 


2mp 


Z = 


Sr^V^pc^ 




IV? 



x\ 



From here on we will refer to p as p for the sake of simplicity. We also define a new 
dependent variable 

p _ dA 
^ dp ' 
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so that we get a system of four first order equations 

A — cos X 

Z 

F 

P 

—Zsmx- (15) 

The equations can also be written as a fourth order equation in the dependent 
variable 

rf' / (d\ . dx\\ ^ d ( (d\ . dx\\ . „ 

Much of the qualitative behaviour of solutions to our system will be determined by the 
behaviour of the dependent variables in the vicinity of the two singular points p = and 

p = CO. 

4 Behaviour Near And oo 

The variable Z is the charge density which we assume to be non-negative (In fact, it is 
straightforward to show that Z > => Z > 0. See proof of Lemma 2.) F is the integral 
of the charge Z per unit ring, that is, F{p) —F{0) is the total charge within a radius p. It 
is reasonable, therefore, to restrict our attention to solutions where F remains bounded 
as p approaches oo. Much of the behaviour of the solution can be characterised by the 
values of F at either end of the domain. This behaviour can best be summarized in the 
following two Lemmas. 

Lemma 1 Suppose {x,F,A,Z) is a solution to Equation(^^ on I = (0,pi), for some 
Pi,0 < pi < 1. Suppose also that Z > is continuous and bounded on /. Then, 

(i) F is on I and has a well defined, finite limit as p ^ 0. Z has a well-defined 
limit as p ^ 0. 

(a) if F{0) 7^ then A is unbounded as p — >■ 0. In particular, A = fi(p) ln(p), where Q 
is and bounded on I, Q ^ F{0) as p 0. Also, x is bounded as p ^ 0. 



(a) 
(b) 
(c) 
id) 



dx 
dp 
dF_ 
dp 
dA 

dp 
dZ_ 
dp 
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Proof 

(i) Let < pi < p2 < 1- 

From (c) , ^ = Z, so F is and \F{p2) - F{pi)\ < M(p2 -pi) where M = supiZ. 
Letting pi, P2 — > shows that F has a well-defined limit. Fixing p2, and letting pi ^ 
shows that this limit is finite. A similar argument using (d); ^ = —Z sinx < M, 
shows that Z has a well-defined limit. 

(ii) Let A = Q\np. A is on / since, from (c), ^ = therefore Q, is C^. We also 
have 

F = p\np^ + n (17) 
dp 



and 



From (0), (since Inp < on /), 



Integrating, 



where 



That is, 



d ( ,n(if2\ 
Mlnp < — p Inp^— < 0. 
dp \ dp J 



M (pi(lnpi - 1) - p(lnp - 1)) < ci - p(lnp)2^ < 0, 

dp 



ci = pi(lnpi)2^(pi). 

dp 

M(pi(lnpi-l)-p(lnp-l)) dn ^ 

< pm p—j— < 



In p dp Inp 

dn 

dp 



Letting p — » 0, we see that p In p^ — > 0. Hence, from (|T^ we see that f2 — > -F(O) 



as p — s> 0. 

Now, since Q is bounded on / put C2 < < C3, for some constants C2 and C3. From 
(a), 

A-l<^<A + l. 
dp 

Thus, 

C2 In p - 1 < ^ < C3 In p + 1. 
dp 

Integrating, we have 

X{p) < x(Pi) + C2(pi(lnpi - 1) -p(lnp- 1)) - (pi -p) 
X{p) > x(Pi) + C3(pi(lnpi - 1) -p(lnp- 1)) + (pi -p). 
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which bounds x{p) as p — > 0. 



□ 



Lemma 2 Suppose {x,F,A,Z) is a solution to Equation on p E (0, oo). Suppose 
also that Z > with F continous and bounded on the interval. Then 

(i) If F{pi) > for some pi G [0, oo), then x oo, A oo, and F ^ oo as p ^ oo. 

(a) If F < on (0, oo) then F ^ as p oo. In addition, if A and Z have well-defined 
limits as p ^ oo then Z ^ and A A^o as p ^ oo, with —1 < A^o < 1. 

Proof 

(i) We firstly show that, given the nondegeneracy condition, u^va a.e., that Z > 
=> Z > for all p E I. Suppose Z{pi) = for some pi G /. Now j-p{ZeP) > 
0, using (d). Integrating from p to pi (p G (0,pi)), shows that Z{p) = for all 
p G (0,pi). A similar argument using ■j^{Ze~^) < shows that Z{p) = for all 
p G (pi, oo). Thus, if Z{pi) = 0, then Z = on I. But this possibility is excluded 
by the non- degeneracy condition. Therefore, Z > on (0, oo). This shows that Z 
does not have isolated zeroes. 

F is (strictly) monotonic increasing everywhere, since from (h) ^ = Z and Z > 0. 
If -F(pi) > for some p2 > pi, F{p2) > 0. So, for p G (p2, oo) 

dA _F ^ F{p2) 
dp p ~ p 

Integrating, 

A{p) > A{p2) + F(p2) In 

and A — > oo as p — >■ oo. 

Also, from (a), ^ > A — 1, which implies, after integration, that x ^ oo as p ^ oo. 

Now on (p2, oo), F is bounded below by F{p2). Let us assume that F is bounded 
above, i.e. < F{p2) < F{p) < (3, a. constant. Consider 

d ( Z{A -cosx) \ _ ZFjl + cosx) . . 

dp[ A+l I p{A + lY ■ ^ ^ 



8 



Now for p2 large enough, A > 1 on {p2, oo) and so 

Z(A-cos x) ^ Z{p2) {A{p2) - cosx(p2)) ^ „ 

'a ; > C4 = — r ; > 0. 

A + l A{p2) + l 

Thus, 

Z(A-cosx) > C4(A + 1). (20) 

Now consider 

V{p) = Z{A-cosx)- — - / -izr^da. 
^ = 0, so V{p) = V{p2), constant on (p2, 00). But 

F{p2f 1 ^, ,2 r da 



Vip) > Z{A - cosx) - - \nP2?f 



2p 2' ^-^^^ Jp, 

and so, using (pO|), since A ^ 00 as p ^ cxd, then V{p) — > cxd as p cxd. Clearly p2 
can be chosen at a point where V^(p2) is finite. We have a contradiction. F cannot 
be bounded by a positive real number. That is, F —>■ 00 as p 00. 



(ii) We assume there exists a constant t] s.t. F{p) < r/ < for all p G (0, 00). Then, 

dA 
dp 

A<A{p2)+vH—), 

P2 



since ^ = F/p, 



and so A ^ —00 as p — >■ 00. Similarly x ~^ ^00 as p 00. 
Now 



d f Z{A - cosx) \ _ ZFi-l + cosx)) (n^ 

dp[ A-1 )~ p{A-iy - 



Therefore 



Z{A - cos x) . _ Z{p2) (A(p2) -cosx(p2)) ^ ^ 
— '^5 — 77 — ^ ; ^ u. 



A-1 - ' A{p2)-1 
for p2 large enough, since ooasp^oo. Therefore 

Z(v4-cosx) < C5{A-1). 
Hence Z{A — cosx) —00 as p — > 00. 

Now, using a similar argument to that in (i), using V{p), we have a contradiction. 
Hence, F cannot be bound above by 77 < 0. But ^ = Z > 0, thus F ^ as 
p — » 00 (otherwise, we have case (i)). The proof of Z — as p — 00 is handled in 
a similar manner. The bounds on A^^ are also easily established using arguments 
of this type. For example, assume Aoo < —1 and use (|19|) and (|2l|) . 

□ 
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5 Numerical Solutions 



5.1 Charged Wire with Two Dense Rings 

As stated in Lemma 2 (ii), if F < on J, then F ^ as p — oo. To obtain a numerical 
solution of this type, we expand x p ^ind solve for the coefficients of Xy ^i F) ^ ^is 
p — s> oo. There is only one such solution, the coefficients in the expansions being uniquely 
determined. 

- TT 3^ 1 1 9247 1 
^ ~ ^~ p^ 6p3 ^ 67320p5 

15 1 35 1 77425 1 
A = H + ... 

2 p2 8 p4 26928 p6 

1 35 1 77425 1 
fp = _i c; I l_ 

p2 2 p4 4488 p6 

1 1 77425 1 , , 

Z = 30^-70^ + —- — . + ... . 22 

p3 p5 748 p7 ^ ' 

Using these expressions, we evaluate our initial values (near infinity). We used a multi- 
step differential equation solver (from the NAG Fortran Library (Mark 16)0), interfaced 
with MATLAB (0)). The numerical solution is given in Figure 1. 

As stated in Lemma 1, x remains finite as p ^ 0. We can calculate the value of x(0) 
numerically: x(0) = —12.17 noting that x ^ tt as p ^ oo. From (d) we have 

dZ 

— = when x = 
dp 

From Figure 1 we see that x is monotonic and bounded between —An and vr. Z has four 
critical points for p G (0, oo). These occur successively as maxima and minima. Z has 
two finite values of p, pi and p2, say, for which the charge density is a local maximum. See 
Figure 1. Numerically, after converting to units, in which p is measured in Ac = ^, the 
reduced Compton wavelength, pi ^ •12Ac, and p2 ~ .46Xc. In the full three-dimensional 
picture, this corresponds to two densely charged rings around the z axis. 

Plotting A against In p, in the inner region (inside the dense rings) we obtain a linear 
plot, in keeping with Lemma l{ii): 

y4^F(0)lnp for p « pi. 

This is just the standard cylindrically symmetric solution to vacuum Maxwell equations - 
the potential due to a charged wire (along the z axis). Note that, F{0) = ^ corresponding 
to a charge per Compton wavelength of 146e. Coupling the Maxwell and Dirac equations, 
then, has the effect of surrounding the charged wire with two dense rings of charge. 
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Charge on ring radius r 



chi vs r 




Figure 1: Numerical solutions to the static D-M equations, with cylindrical symmetry 
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Charge on ring radius r 




Figure 2: Numerical solution to the Dirac equation aound a charged wire, treated as an 
external field 

5.2 Comparison with the Linearized Theory 

In this section we consider, by way of comparison, the solution to the Dirac equation in 
an external field generated by a charged wire. We now use the decoupled equations, i.e. 
the Dirac equation with the vacuum Maxwell equation. 

In the static cylindrically symmetric case, the solution to the vacuum Maxwell equation 
is A = Co In p. Using the same notation as in the previous section, we choose cq = -^(0) < 0. 
We can now solve the Dirac equation, using this (external) potential. Our equations are 



(&) 



dp 
dZ 
dp 



— 10.79 Inp — cosx 
— Zsinx. 



(23) 



For comparison purposes, we solve these linearized equations with the same initial condi- 
tions used for our earlier solution. The numerical solution was obtained using the same 
methods as in the previous section and are shown in Figure 2. 

For small values of p the current density, Z, behaves in a similar manner (to the 
solution of the non-linearised D-M equations), having the first peak at pi ^ .llAc- As 
p ^ oo the (logarithmic) potential is unbounded, in contrast to the nonlinearised case, 
in which — 1 < A^o < 1 (see Lemma 2). It is clear from equation (|2^(a)) that x ^ — oo 
as p — ^ oo. As such, (from ^(6)), there will be an infinite number of oscillations in the 
charge density variable Z, which, however, remains bounded. We have lost the localization 
that was apparent in the full D-M equations. 
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Figure 3: Charge distribution of the Dirac-Maxwell field with cyhndrical symmetry evi- 
dencing a highly localised charge density around a central charged wire. Total charge (per 
unit length) is finite. 
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Figure 4: Charge distribution of the Dirac equation surrounding a charged wire, treated 
as an external field. No localisation of the charge distribution is apparent. The total 
charge (per unit length) is unbounded. 
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5.3 Unbounded Solutions 



As stated in Lemma if F(0) 7^ then A is unbounded as p — > 0. We can also 

look for solutions where F(0) = 0. Lemma 2 then tells us that x, ^4, F — > 00 as p — > 00. 
Since F{p) is the charge within a radius p, these solutions are of less interest, as the total 
charge, F{oo) is necessarily unbounded. 

To find a numerical solution of this type, we expand x iiito a Taylor Series in p and 
solve for the coefficients. 

As an example of this type of solution we set A{0) = 1 and Xo = 0. The first few 
terms of these and the remaining variables are given below. 

_ J_ J_ 2 1 3 1 4 
^ ~ 100^^100^ ^ 60000^ ^ 45000 

_ 101 J_ ]_ 3 1 4 

100^50^ 90000^ 240000^ 
F = —p p* + ... 

50'^ soooo'^ eoooo'^ 

Z = p^ p^H p^ + ... . (24) 

50 10000^ 15000^ 4000000000^ ^ ^ 

In Figure 5 we show the solution which has these initial values. 



6 Conclusions 

An examination of the behaviour of the solutions to the static Dirac-Maxwell Equations 
with cylindrical symmetry shows that, in one class of solution, we have a highly localised 
charge density around a central charged "wire" . The total charge (per unit length of the 
"wire" ) is finite. 

When the equations are decoupled (by ignoring the effect of the "self-field" upon the 
Maxwell field) all localisation is lost. The total charge (per unit length) of the Dirac field 
is, in this case, unbounded. 

Although the cyhndrical case is of hmited interest physically, the resulting o.d.e.s are 
amenable to a detailed descriptive analysis which corrorborates the numerical results. 

The same localisation due to the inclusion of a nonlinear coupling, will be more in- 
teresting in the axially symmetric case, an analysis of which will be in a forthcoming 
publication. 
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